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“EXPLORING  NOVEL  FORMS  OF  SUPERCONDUCTIVITY....” 


1 .  New  Growth  mode 

First  we  focused  on  developing  a  growth  mode  that  is  capable  of  producing  a  single  unit  cell 
thin  superlattices  of  LaNi03/LaAL03  (LNO/LAO  thereafter).  High-quality  epitaxial  LNO  ultra- 
thin  films  (10-15uc)  were  grown  on  STO  (001)  and  LAO  (001)  single  crystal  substrates  by 
pulsed  laser  deposition  (PLD)  with  in-situ  monitoring  by  Reflection  High  Energy  Electron 
Diffraction  (RHEED).  Due  to  the  lattice  mismatch,  NNO  (3 .80 A)  is  subjected  to  a  tensile  strain 
of  2.6  %  on  STO  (3 .90 A),  whereas  a  compressive  strain  of  -0.3  %  is  applied  on  LAO  (3.79A). 
The  used  LAO  substrates  had  mixed  terminations  as  received,  while  the  STO  substrates  were 
prepared  by  our  recently  developed  chemical  wet-etch  procedure  (Arkansas  treatment)  to  attain 


Time  (s) 


Figure:  (a)  A  representative  RHEED  specular  intensity  evolution  during  the  continuous  slow  growth  and  the  interrupted 
growth  with  fast  laser  ablation.  Inset:  the  intensity  oscillation  during  the  ”laser-on”  period  (highlighted)  of  the 
interrupted  growth.  AFM  scans  on  the  surface  of  a  Ti02-  terminated  SrTi03  substrate  before  the  growth  (b),  after  the 
continuous  slow  growth  (c)  and  the  interrupted  growth  (d). 


an  atomically  flat  Ti02-terminated  surface  with  a  minimal  number  of  surface  bound  electronic 
defects.  The  growth  temperature  was  in  the  range  of  670-700  C,  while  the  oxygen  partial 
pressure  was  maintained  at  100-120  mTorr.  After  deposition,  to  maintain  the  proper  oxygen 
content,  samples  were  post  annealed  for  30  minutes  and  cooled  down  to  room  temperature  in  1 
atmosphere  of  ultra-pure  oxygen. 
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To  maintain  the  perfect  epitaxy,  LBL  growth  is  achieved  by  the  interrupted  deposition, 
which  requires  a  rapid  laser  ablation  cycle  (up  to  30Hz)  followed  by  a  prolong  time-delay 
between  two  successive  unit  cells.  Figure  1(a)  and  Fig.  2  shows  a  representative  time-dependent 
RHEED  specular  intensity  (RSI),  where  the  full  recovery  of  RSI,  characteristic  of  the  perfect 
LBL  growth,  can  be  seen  after  each  unit-cell  layer.  Smooth  surface  morphology  is  also  evident 
from  the  well-defined  spots  of  the  (00)  specular  and  (01)  and  (0-1)  off-specular  reflections  with 
streaks  in  the  resulting  RHEED  pattern.  AFM  imaging  showed  that  the  sample  surface  is 
atomically  flat  with  preserved  vicinal  steps  and  typical  surface  roughness  <  80  pm. 
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Figure  2.  Surface  morphology  and  structural  characterization  of  the  LNO  films,  (a)  Time-dependent  RHEED 
specular  intensity  (RSI),  where  the  full  recovery  of  RSI,  characteristic  of  the  perfect  layer-by-layer  (LBL)  growth,  is 
observed  after  each  unit-cell  layer  [indicated  by  the  shaded  (grey)  area].  Smooth  surface  morphology  is  evident  from 
the  well-defined  spots  of  the  specular  and  off-specular  reflections  with  a  characteristic  to  quality  surface  streaks 
pattern  shown  in  inset,  (b)  AFM  imaging  confirms  that  the  sample  surface  is  atomically  flat  with  well-preserved 
vicinal  steps  and  typical  surface  roughness  <80  pm.  (c)  A  reciprocal  space  map  around  (222)  reflection  of  10  u.c. 
LNO  film  grown  on  STO.  A  position  of  the  reflection  is  in  perfect  registry  with  the  substrate.  The  absence  of  diffuse 
intensity  around  the  peak  testifies  for  the  full  epitaxy,  (d)  High  resolution  TEM  image  at  the  film-substrate  interface 
(shaded)  also  support  excellent  sample  epitaxy. 

To  elaborate  the  role  of  the  laser  frequency  f  during  the  growth  of  each  unit-cell  layer,  Fig.  3 
shows  a  schematic  of  the  evolution  of  the  growth  dynamics.  The  nucleation  and  crystallization  in 
pulsed  laser  deposition  (PLD)  are  primarily  governed  by  two  time  scales;  one  is  the  length  of  the 
”flux-on”  period  tF  during  which  the  dense  plasma  plume  typically  enforces  high  supersaturation 
and  nucleation  of  small  islands  on  the  order  of  one  structural  unit,  while  the  other  is  the 
’’relaxation”  period  tR  where  the  absorbed  atoms  and  small  islands  reach  stabler  state  by 
diffusing  on  the  surface  to  form  larger  islands  or  migrate  to  the  edges  of  larger  islands.  tF  is 

relatively  short  on  the  scale  of  lO-^  to  10-^  s.  In  contrast,  tR  is  about  10_~  to  10-^  s  and  much 
longer.  When  tR  «  1/  f  ,  the  small  islands  always  migrate  and  reach  a  steady  state  before  the 
material  is  delivered  again  by  the  next  pulse  (see  Fig.  3(a)).  In  this  case,  the  LBL  growth  relies 
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on  the  condition  that  the  absorbed  atoms  have  enough  mobility  to  find  and  fill  the  voids  of  the 
incomplete  growing  layer.  Otherwise,  the  growth  will  transform  into  a  3D  island  fashion  due  to 
stacking  of  several  incomplete  layers  at  the  growth  front.  In  contrast,  when  tR  ~  1/  f  ,  due  to  the 
transitional  kinetics  imposed  by  the  rapid  laser  pulsing,  the  concentration  of  the  small  nucleating 
islands  does  not  fully  decay  between  two  consecutive  pulses  (see  Fig.  3(b)).  In  this  case,  the 
growing  layer  is  formed  by  homogeneous  coalescence  of  small  islands  instead  of  isolated  large 
islands,  reflected  in  the  better  recovery  of  the  RHEED  specular  intensity.  If  one  is  able  to 
decrease  1/  f  to  the  order  of  tF,  the  system  will  see  only  an  average  flux  of  vapor,  and  the  growth 
will  fall  into  a  continuous  MBE-  like  regime  (see  Fig.  3(c)).  Finally,  since  there  are  fewer  and 
fewer  voids  as  the  layer  coverage  improves,  some  atoms  may  remain  on  the  surface.  In  this  case, 
an  interruption  on  the  laser  ablation  can  be  introduced  to  delay  the  delivery  of  the  next  layer  so 
as  to  wait  for  the  layer  completion.  This  corresponds  to  the  slow  recovery  of  the  RHEED 
specular  intensity  in  the  ”laser-off  ’  period. 


Figure  3:  Timing  diagram  showing  the  modulated  flux  and  relaxation  of  the  island  density.  During  the  low- 
frequency  pulsed  growth  (a),  the  physical  vapor  is  delivered  to  the  substrate  surface  where  it  is  relaxed.  After 
increasing  the  pulse  rate  (b),  the  re-  laxation  process  becomes  more  and  more  delayed.  As  the  frequency  (c)  reaches 
the  quasi-constant  state  with  average  flux,  the  large  degree  of  supersaturation  overcomes  the  nucleation  barrier 
result-  ing  in  the  rapid  nucleation  from  the  supercritical  vapor. 

2.  Effect  of  Polar  Mismatch 

Recently,  the  problem  of  polar  interfaces  has  been  drawing  enormous  attention  for  building 
interface-  controlled  devices,  since  a  polar  crystal  structure  grown  on  a  non-polar  substrate 
would  be  unstable  due  to  the  depolarizing  fields.  An  example  is  showed  in  Fig.  4(a)  where  two 

perovskites  (e.g.  LAO  and  STO)  with  different  charge  structures,  such  as  a2+B^+02  and 

A2+b4+o2'  ,  are  interfaced  together.  The  +1 — 1 — hi — 1  alternating  stacking  sequence  on  the 
LAO  side  would  resulting  in  a  diverging  electrostatic  potential  or  a  diverging  electric  dipole 
moment,  as  the  LAO  slab  becomes  thicker  and  thicker.  To  avoid  this  unstable  divergence,  some 
mechanism  must  set  in  to  compensate  the  potential  jump  by  forcing  electronic,  structural  or 
chemical  (non-stoichiometry)  changes.  In  particular,  the  active  work  on  this  LAO/STO  system 
has  suggested  an  electron  transfer  to  the  interfacial  Ti02  that  decreases  Ti  valence  by  0.5  at  the 
expense  of  0.5e  charge  in  the  terminating  atomic  plane  (see  Fig.  4(b)).  This  mechanism  has  been 
intensively  discussed  to  account  for  the  observed  2D  electron  gas  at  the  LAO/STO  interface, 
although  additional  effects  such  as  alloying  and  ionic  displacement  have  been  suggested  at  the 
interface. 
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a)  b)  c) 


Figure  4:  Cartoon  pictures  of  polar  interfaces  of  LAO/STO  before  charge  transfer  (a),  after  charge  transfer  (b),  and 
with  some  of  the  LAO  layers  replaced  by  LNO  in  a  SL  structure  (c).  V  denotes  the  electron  potential. 

A  simple  fundamental  question  one  can  ask  is  that  what  if  some  of  the  ions  at  the  LAO  side  are 
replaced  by  elements  that  have  multiple  valence  states.  More  interestingly,  to  test  the  charge- 
transfer  mechanism,  we  may  replace  some  the  LAO  layers  with  LNO  (see  Fig.  4(c))  where  Ni  is 
then  supposedly  converted  from  3+  toward  4+  valence  state  by  giving  up  electrons.  In  this 

section,  we  will  show  that,  instead  of  Ni^+,  an  unusual  valence  state  of  Ni~+  is  observed  in  ultra- 
thin  LNO/LAO  SLs  grown  on  non-polar  (or  weakly  polar)  (001)  Ti02-terminated  STO 
substrates.  This  is  in  marked  contrast  to  identical  SLs  grown  on  polar  LAO,  which  always 
exhibit  the  3+  valence  of  the  bulk-like  LNO.  To  clarify  the  origin  of  these  phenomena,  we  have 
tracked  the  changes  of  the  Ni  valence  as  a  function  of  the  SL  thickness  and  developed  an 
experimental  procedure  to  monitor  and  control  the  electronic  state  of  a  transition  metal  ion. 

3.  NOVEL  ELECTRONIC  PROPERTIES 


To  get  insight  into  the  electronic  structure,  we  carried  out  detailed  soft  x-ray  absorption 
spectroscopic  measurement,  which  is  particularly  useful  in  probing  layers  buried  within  SLs  due 
to  the  local  and  element-specific  characters.  Spectra  were  acquired  in  the  bulk-sensitive 
fluorescence  yield  mode  at  the  Ni  L-edge  at  the  4ID-C  beamline  of  the  Advanced  Photon  Source. 
The  normalized  spectra  taken  at  300  K  are  shown  in  Fig.  5.  The  main  peaks’  energies  and  the  line 
shape  are  similar  to  the  bulk  data  of  ReNi03  ,  indicative  of  a  Ni3+  valence  state  with  the  3d7 
low-spin  configuration.  On  the  other  hand,  as  n  decreases,  the  splitting  of  the  double  peaks 
(highlighted  in  green)  at  the  L3-edge  becomes  more  pronounced,  especially  from  n  =  10  to  5. 
The  corresponding  trend  at  the  L2-edge  emerges  as  a  weak  shoulder  at  about  870.5  eV.  Notice 
that  this  multiplet  splitting  is  absent  in  bulk  LNO,  but  reminiscent  of  the  evolution  as  entering 
the  charge-ordered  insulating  state  in  bulk  ReNi03.  The  emergence  of  this  multiplet  splitting  in 
the  SLs  even  with  a  relatively  thick  LNO  layer  (n  =  10)  denotes  enhanced  localization  and  hints 
at  a  MIT  and  a  latent  tendency  to  charge  ordering. 
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X-ray  linear  dichroism  (XLD)  (iEiab  -  Ieic)  was  then  measured  to  acquire  the  eg  orbital 
occupation.  Spectral  weight  transfer  in  XLD  shown  in  Fig  5  between  the  main  peak  and  the  low- 
energy  multiplet  is  indicated  by  the  opposite  signs,  implying  the  Oh  symmetry  of  the  charge 
distribution  is  reduced  likely  to  D4h.  The  main  peak’s  negative  XLD  signifies  more  holes 
populate  the  dz2  orbital,  resulting  from  tensile  strain  and  interfacing  with  LAO.  As  n  increases, 
the  XLD  is  rapidly  suppressed,  and  becomes  practically  zero  for  n  =  10.  The  origin  for  this  can 
be  twofold.  First,  at  large  thickness,  the  inherent  orbital-lattice  interaction  of  LNO  tends  to 
preserve  the  eg  degeneracy  in  response  to  tensile  strain.  Second,  in  a  localized-itinerant 
crossover,  such  as  Mott  transition,  the  kinetic  energy  gain  favors  putting  electrons  in  orbitals  that 
participate  in  band  formation  as  in  manganites.  In  our  case,  that  is  the  dz2  orbital  with  increasing 
n. 


Figure  5.  X-ray  absorption  spectra  of  (LNO)n/(LAO)3  SLs  at  the  L3  (left  panel)  and  L2  (right)  edges  for  n  =  3,  5, 
and  10  (top).  The  bottom  shows  the  theoretical  calculations  done  for  energy  (eV)  Ni06  cluster  for  a  charge-transfer 
energy  of  0.8  (red)  and  1.6  (blue)  eV.  The  shaded  (green)  areas  indicate  where  the  most  significant  changes  occur 
due  to  the  changes  in  covalency. 


To  identify  the  ground  states,  DC  transport  measurements  from  400  K  to  2  K  were 
performed  in  a  commercial  physical  properties  measurement  system.  The  temperature  dependent 
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conductivity,  o(T),  of  the  SLs  is  shown  in  Fig.  6.  A  transition  from  metallic  to  insulating 
behaviors  can  be  clearly  seen  with  decreasing  n.  The  n-dependence  of  a  (the  left  inset  of  Fig.  6) 
also  shows  orders  of  magnitude  change  across  the  series.  The  n  =  5  and  10  SLs  are  metallic  at 
high  temperature  and  reach  conductivity  maxima,  a  max?  Tmax—  120  K  and  10  K,  respectively. 
For  the  n  =  10  SL,  Ao(T )  (o(Tmax)  -  o(2  K))  is  only  ~  52.5  S/cm,  relatively  small  compared  with 
o(2  K)  (~  1%).  This  do/dT  >  0  region  is  attributed  to  quantum  correction  for  metal.  In  contrast, 
Ao(T  )  (-381 .6  S/cm)  for  n  =  5  is  similar  to  o(2  K)  and  almost  twice  of  the  increase  in  o(T)  from 
400  K  to  Tmax.  This  behavior  implies  the  system  is  on  the  verge  of  MIT. 

To  quantify  the  transport  properties,  we  fitted  the  data  to  various  temperature  dependent 
behaviors.  The  low  temperature  part  of  the  n  =  3  SL  is  found  to  be  well  fitted  with  the  2D 
variable  range  hopping  behavior  (lno  °c  -T~1/3)  up  to  20  K.  It  switches  to  3D  variable  range 
hopping  (lno  -T~1/4)  up  to  around  80  K.  The  presence  of  hopping  conductivity  reveals  that  the 
system  is  on  the  insulating  side  of  the  MIT  when  n  =  3.  On  the  metallic  side,  the  resistivity  of  the 
n  =  10  SL,  q,  follows  a  T  -linear  behavior  above  200  K  and  then  enters  a  crossover  region  before 
switching  to  a  T1 5-dependence  below  110  K.  The  low  temperature  T1 5-dependence  may  be 
related  to  proximity  to  an  itinerant  to  localized  transition.  Below  Tmax,  while  quantum 
correction  often  consists  of  multiple  contributions,  fitting  with  a  complete  expression  is 
inconclusive  due  to  the  limited  range  of  data.  Instead,  fitting  with  a  simple  power-law  relation  (a 
oc  Tm)  yields  m  ~  1/3  for  data  below  5.5  K.  It  is  interesting  to  note  that,  while  m  =  1/2  arising 
from  electron-electron  interaction  can  describe  quantum  correction  in  many  correlated  metals,  m 
is  found  to  decrease  from  1/2  to  1/3  in  a  number  of  metallic  perovskite  oxides,  including  LNO, 
when  approaching  the  MIT  from  the  metallic  side  upon  reducing  the  carrier  density.  For  the  n  =  5 
SL,  q  is  also  linear  with  temperature  above  270  K.  The  region  between  that  and  Tmax  appears  to 
be  a  crossover  as  no  simple  power-law  dependence  is  found. 

Unlike  the  n  =  3  and  10  SLs,  the  do/dT  >  0  region  of  the  n  =  5  SL  cannot  be  accounted  by 
hopping  conductivity  or  quantum  correction.  Imposing  the  Tm  power-law  relation  gives  rise  to  a 
strongly  temperature-dependent  m,  which  increases  continuously  to  about  0.9.  This  interesting 
behavior  is  probably  related  to  the  fact  that  Ao  ~  Oo,  indicating  that  all  the  interactions  involved 


6 


in  the  dimension-control,  such  as  correlation  effects  and  the  bandwidth,  are  of  a  similar  strength. 
The  n  =  5  SL,  thus,  falls  into  a  critical  region  of  the  MIT. 


Figure  6:  (a)  Conductivity  of  (LNO)w/(LAO)3  SLs  versus  temperature.  Inset:  conductivity  versus  n  at  10  K.  The 

dashed  line  is  a  guide  for  eye.  (b)  Resistivity  as  a  function  of  T  ^  ^  at  n  =  10.  The  red  solid  line  is  a  guid  for  eye.  (c) 
1/4 

dlna  as  a  function  of  1  IT  for  the  n  =  3  SL  shows  the  low  temperature  three  dimensional  variable  range  hopping 
behavior. 

4.  Results  and  Discussion  on  electronic  properties 

High-quality  epitaxial  [(LNO)n/(LAO)3]/v  SLs  (n  =  3,  5  and  10  with  N  =  10,  4  and  2 
respectively)  were  grown  by  laser  molecular  beam  epitaxy  (MBE)  as  described  above.  +2.2% 
tensile  strain  is  applied  with  Ti02-terminated  (001)  SrTiC>3  (STO)  substrates.  To  explicate  the 
evolution  of  the  ground  state  due  to  quantum  confinement,  DC  transport  measurements  from  400 
K  to  2  K  were  performed  by  using  the  van  der  Pauw  method.  The  obtained  temperature- 
dependent  conductivity  o(T )  of  the  SL  series  is  shown  in  Fig.  6(a).  A  transition  from  metallic  to 
insulating  behavior  can  be  clearly  seen  with  decreasing  n.  The  inset  of  Fig.  6(a)  also 
demonstrates  orders  of  magnitude  change  in  a  across  the  series  taken  at  10  K.  As  seen,  in 
analogy  to  the  bulk  LNO,  the  n  =  10  SL  exhibits  a  highly  metallic  behavior  ascribed  to  the  three 

3/2 

dimensional  regime.  On  the  other  hand,  unlike  the  bulk,  which  is  known  to  follow  a  T 

dependence  at  room  temperature  and  switch  to  a  low-temperature  T  Fermi  liquid  behavior,  the 

resistivity  Q  at  n  =  10  switches  from  a  T-linear  behavior  above  200  K  to  a  T  power  law  (see 
Fig.  6(b)).  Note  that,  while  a  T- linear  dependence  was  also  found  in  the  high-temperature 
metallic  phase  of  other  members  of  the  bulk  RNi03  family  who  have  a  charge-ordered  insulating 

ground  state,  the  low -temperature  T  behavior  of  Fermi  liquid  electrons  has  been  attributed  to 
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scattering  by  bond-length  fluctuations  due  to  proximity  to  a  Mott-type  MI  transition  from  the 
itinerant  side.  These  phenomena  indicate  enhanced  electronic  correlations  in  the  quantum 
confinement  regime  and  stronger  vibronic  fluctuations  associated  with  charge  ordering. 

The  intermediate  thickness  n  =  5  SL  still  shows  a  clear  metallic  behavior  at  high  temperatures 
but  reaches  a  conductivity  maximum  at  Tmax=  120  K.  While  the  metallicity  of  LNO  films  under 
large  tensile  strain  is  known  to  cease  when  the  thickness  is  <  3nm  ,  it  is  remarkable  that  the  high- 
temperature  metallic  behavior  is  well  maintained  in  such  a  thin  slab  (<  2nm).  This  result 
unambiguously  highlights  the  exquisite  role  of  the  confinement  structure  with  interfaces  in 
defining  the  transport  properties  of  SLs  as  compared  to  the  bulk  LNO  and  single  layer  films. 
Notice,  the  conductivity  a  here  is  also  well  below  the  Ioffe-Regel  limit  (~  2000  S/cm), 
manifesting  the  modification  of  the  correlation  nature  of  the  observed  behaviors  and  the  minor 
role  of  dis-  order.  Below  Tmax,  o(7max)  -  o(2K)  (-381.6  S/cm)  is  as  large  as  o(2K)  and  almost 
twice  of  a(rmax)  -  o(400  K).  Neither  hopping  conductivity  nor  quantum  corrections  is  found  to 
be  able  to  account  for  this  daldT>  0  region.  This  transitional  behavior  explicitly  indicates  that  all 
the  involved  interactions  and  relevant  energy  scales,  like  electron  correlations  and  the  bandwidth, 
are  of  a  similar  strength  and  strongly  competing  with  each  other.  Thus,  n  =  5  falls  into  a  critical 
region  of  the  confinement-controlled  MIT.  Upon  further  decreasing  n  to  3,  do/dT  >  0  at  all 
temperatures  and  a  variable  range  hopping  behavior  was  found  below  100  K,  indicating  that 
finite  density  of  states  remains  within  the  opening  gap  during  the  dimensional  crossover. 
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Figure  7:  Layer  CDMFT  results  for  the  multilayer  Hubbard  model  with  various  thickness  n.  The  transfer  anisotropy 
6  is  set  at  0.2.  (a)  Average  spectral  function  for  U  It  -  6.  (b)  Sheet  conductance  a  versus  n  for  Ult  -5.5,6  and  6.5. 


To  elucidate  the  origin  of  the  MI  transition  induced  by  quantum  confinement,  we  compare  the 
dimensional  crossover  with  the  layer  extension  of  cellular  dynamical-mean-field  theory  (layer 
CDMFT)  which  takes  into  account  the  competition  between  the  short-range  correlations  and  the 
reduced  bandwidth.  To  make  the  problem  tractable  and  yet  relevant  for  the  SLs,  we  considered 
the  multi-layered  single-band  Hubbard  model  for  Ni  3d  and  simulated  the  high-energy  A1  3s 
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states  as  vacuum.  Note  that  CDMFT  was  recently  applied  to  explain  the  opening  of  a  correlation- 
induced  pseudogap  in  SrV03  thin  films.  The  effective  impurity  model  for  each  layer  consisting 
of  2x2  plaquette  coupled  to  8  bath  orbitals  and  is  solved  by  using  the  Lanczos  exact 
diagonalization  method.  The  layer  CDMFT  employed  here  goes  beyond  that  previous  study  by 
including  short-range  correlations  within  each  monolayer  to  capture  the  physics  of  pseudogap. 
The  short-range  correlation  is  antiferromagnetic  as  recently  suggested  in  a  similar  system  with 
the  planar  dx2_y2  orbital  configuration.  To  neutralize  the  CDMFT’s  overemphasis  on  the  in¬ 
plane  correlations,  we  introduce  a  small  anisotropy  between  the  in-plane  tjn  and  out-of-plane  tout 
hopping  integrals  as  tin  =  t  ( 1  -  8)  and  tout  =  t  (1  +  28)  so  that  tout  competes  with  the  in-plane 
correlations.  The  bandwidth  of  the  bulk  non-interacting  system  is  unchanged,  W  =  12 1. 

Since  bulk  LNO  is  a  metal  we  consider  moderate  values  of  U/t  (U  =  6 1).  Fig.  7(a)  presents  the  n- 
dependent  evolution  of  the  spectral  function.  As  clearly  seen,  the  correlation  gap  in  the  two 
dimensional  limit  is  filled  up  and  becomes  a  pseudo-gap  with  increasing  n,  lending  theoretical 
support  to  the  notion  that  the  experimentally  observed  MI  transition  arises  from  the  dimensional 
evolution  of  a  quantum  confined  correlated  carriers  transforming  from  an  insulator  to  a  metal 
through  a  ’’bad”  metal  phase.  Notice,  the  sharp  quasiparticle  peak  is  absent  at  to  =  0  even  in  the 
thickest  slab  because  out-of-plane  correlations  are  treated  on  a  mean-field  level.  To  see  the 
relation  between  the  (pseudo)  gap  opening  and  the  anomalous  transport  properties,  sheet 
conductance  a  is  computed  by  the  standard  Kubo  formula.  An  imaginary  part  r)  =  0.1  f  was 
introduced  to  the  real  frequency  to  to  represent  the  finite  impurity  scattering.  As  shown  in  Fig. 
7(b),  with  increasing  n,  o  increases  continuously  and  becomes  practically  unchanged  above  a 
critical  thickness  nc  below  which  the  pseudo-gap  is  pronounced.  These  results  reproduce  well  the 
critical  regime  at  nc  ~  5  in  the  transport  measurements  and  strongly  suggest  the  confinement- 
controlled  short-range  electron  correlations  as  the  origin  of  the  observed  MI  transition.  From  the 
comparison  above,  we  can  see  that  the  essential  physics  of  metal-insulator  crossover  by  the 
quantum  confinement  is  captured  by  the  CDMFT.  On  the  other  hand,  as  discussed  earlier,  it  is 
known  that  rare-earth  nickelates  also  have  an  intrinsic  propensity  towards  the  local  charge  dis¬ 
proportionation  upon  entering  the  insulating  phase.  The  result  is  a  complicated  3D  ordering  of 
both  charge  and  spin  which  might  emerge  as  a  competing  ordering.  To  gain  microscopic  insight 
into  the  local  electronic  structure  and  the  effect  of  confinement  on  the  MI  transition,  we  carried 
out  detailed  x-ray  absorption  spectroscopy  (XAS)  measurements.  All  spectra  were  acquired  in 
the  bulk- sensitive  fluorescence-yield  mode  at  the  Ni  L-edge.  The  normalized  spectra  obtained  at 
300  K  (not  shown).  The  La  A/4-edge  which  partially  overlaps  with  the  Ni  L3-edge  was 
subtracted  with  a  Lorentzian  profile.  For  direct  comparison,  the  L-edge  spectra  of  bulk  LNO  and 
SmNi03  are  included.  As  seen,  the  energy  position  and  the  lineshape  of  the  SLs  are  close  to 

those  of  the  bulk  RNi03,  indicative  of  a  Nf®  charge  state  of  the  low-spin  configuration.  On  the 
other  hand,  while  the  multiplet  effects  are  smeared  in  the  metallic  bulk  LNO,  the  L3-edge  spectra 
of  the  SL  show  a  pronounced  two-peak  structure  (shaded  in  green),  closely  resembling  the 
absorption  spectrum  of  the  bulk  SmNi03  ,  which  in  the  bulk  is  a  well  studied  case  for  the 

charge-ordered  (2Ni^-*Ni^+Ni^  )  insulating  state.  As  n  decreases  this  multiplet  effect 
becomes  more  pronounced,  particularly  in  moving  from  n  =  10  to  5.  The  energy  separation  of  the 
two  split  peaks  is  also  clearly  increasing.  The  corresponding  trend  at  the  L2-edge  emerges  as  a 
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developing  low-energy  shoulder  at  about  870.5  eV.  The  observed  evolution  is  reminiscent  of  that 
reported  for  the  bulk  RN1O3  when  crossing  the  MI  boundary,  unambiguously  revealing  the 
carrier  localization  in  the  presence  of  a  developing  correlated  charge  gap  in  the  quantum- 
confined  LNO  slabs  and  a  latent  tendency  to  CO  during  the  3D-to-2D  crossover.  In  accordance 
with  the  observations  from  DC  transport,  the  emergence  of  the  multiplet  even  in  the  relatively 
thick  LNO  slab  (n  =  10)  lends  solid  support  to  strongly  confinement-enhanced  correlations. 
Meanwhile,  when  subject  to  this  enhancement,  the  inclination  to  CO  corroborates  its  role  as  a 
primary  competing  ground  state  not  only  in  the  bulk  but  also  in  the  predicted  heterostructured 
LNO. 

4.  EFFECTS  OF  HETEROEPITAXIAL  STRAIN 


Heteroepitaxial  strain  is  one  of  the  most  important  and  most  extensively  utilized  control 
parameters  in  epitaxy  engineering.  It  is  also  an  inevitable  effect  because  lattice  mismatch  always 
exist  between  two  different  materials.  As  directly  acting  on  the  lattice  of  a  thin  layer  of  material 
by  interfacial  bonding,  epitaxial  strain  is  often  more  efficient  than  elastic  mechanical  stress  and 
the  resulting  property  modifications  can  be  dramatic.  Therefore,  to  rationally  design  a 
heterostructure,  the  lattice  mismatch  is  an  essential  starting  point.  In  this  section,  we  will  discuss 
how  the  bi- axial  character  of  heteroepitaxial  strain  affects  the  orbital  degree  of  freedom  of  d- 
electrons  in  LaNiC>3  (LNO)  ultrathin  films,  since  orbital-lattice  interaction,  e.g.  the  crystal  field 
effect  and  the  JT  effect,  is  critical  to  the  local  electronic  configuration  and  related  long-range 
orderings.  The  second  part  will  go  further  to  show  the  control  of  heteroepitaxial  strain  on  the  p  - 
d  covalency  via  modulating  both  the  self-doping  and  hybridization.  As  a  result,  the  metal- 
insulator  (MI)  transition  and  the  collective  charge-  and  spin-ordered  ground  state  of  NdNi03 
(NNO)  ultra-  thin  films  can  be  fine  tuned  within  a  larger  phase  space. 

5.  Orbital  Response  to  Lattice  Mismatch 

Heteroepitaxy  is  a  powerful  venue  to  modify  orbital-lattice  interactions  in  correlated  materials 
with  strong  electron-electron  interactions  derived  from  transition  metals  with  open  d-shell 
configurations.  However,  efforts  to  rationally  control  properties  that  are  exceedingly  sensitive  to 
small  perturbations  through  the  strong  orbital-lattice  coupling  are  impeded  by  a  poor 
understanding  of  how  the  epitaxial  interfaces  impose  constraints  on  the  orbital  response.  Thus, 
despite  the  recent  progress  in  strain-induced  orbital  engineering,  a  crucial  fundamental  question 
remains:  when  a  single  electron  occupies  a  doubly  degenerate  manifold  of  ci-orbital  in  an 

octahedral  crystal  field,  such  as  in  perovskites  with  Cu~+,  Mir1+  or  low-spin  Ni^+  cations,  how 
do  the  interface-imposed  epitaxial  constraints  select  the  correlated  orbital  response  and  electronic 
configurations? 

The  exceptional  strain-control  of  the  frontier  atomic  orbitals  relies  on  the  susceptibility  of  the 
orbital  occupations  and  their  energy  splitting  to  bi-axial  strain-induced  lattice  deformations.  For 
this  reason,  the  conventional  orbital  engineering  approach  in  perovskite  oxides  is  often 
rationalized  as  follows:  heteroepitaxial  strain  imposes  a  tetragonal  distortion  on  the  film’s 
primitive  unit  cell,  which  results  in  a  modification  of  the  M-0  bond  lengths  of  the  M06 
octahedral  units.  The  bond  distortions  in  turn  reduced  the  crystal  field  symmetry  and  remove  the 
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two-fold  eg  orbital  degeneracy.  Since  the  x  -  y  and  3 z  —  r  orbitals  represent  in-plane  and 
out-of-plane  electron  densities,  respectively,  it  is  generally  anticipated  that  tensile  and 
compressive  strains  should  symmetrically  split  their  energy  levels  but  with  opposite  signs.  The 
bands  derived  from  these  two  orbitals  are  either  lowered  or  raised  relative  to  their  strain-free 
band  center  of  mass,  leading  to  an  orbital  polarization  with  a  finite  filling.  This  symmetric  strain- 
induced  orbital  polarization  (SIOP)  concept  is  routinely  used  to  rationalize  the  orbital  responses 
of  many  complex  oxide  systems,  and  also  plays  a  key  role  in  several  recent  theoretical  proposals 
for  replicating  high-Tc  cuprate  physics  in  perovskite  nickelate  heterostructures.  However,  the 

strategy  for  tuning  the  correlated  ground  states  in  ultra  thin  films  devised  from  this  SIOP  concept 

2  2 

is  violated  more  often  than  it  applies:  for  instance,  doped  manganite  thin  films  exhibit  3 z  -  r 
orbital  polarization  plainly  inconsistent  with  the  SIOP  model  regardless  of  whether  the  applied 
bi- axial  strain  is  tensile  or  compressive.  Ultra  thin  cuprate  bilayers  also  show  variable  critical 
temperatures  depending  on  the  layering  sequence  correlated  with  the  sign  of  the  interface  lattice 
misfit,  and  an  ferromagnetic  order  can  be  switched  on  (or  off)  depending  on  the  strain  state  in 
otherwise  diamagnetic  cobaltite  films.  In  this  section,  we  will  provide  insight  into  the  violation  of 
the  SIOP  model  via  our  investigation  on  the  orbital-lattice  interactions  in  10-u.c.  thick  films  of 

orbitally  degenerate  (e^)  LNO. 

6.  Results  and  Discussion  on  effects  of  strain 

Here  we  focus  on  the  evolution  of  the  orbital  polarization  under  two  representative  strain  states: 
compressive  and  tensile  strain  on  LaA103  (LAO)  and  SrTi03  (STO)  substrates  respectively.  We 
synthesized  high-quality  (OOl)-oriented  10  unit-cell  (u.c.)  epitaxial  LNO  films  by  laser  molecular 
beam  epitaxy  (MBE)  as  described  above.  Single  crystal  LAO  and  STO  were  used  to  apply 
opposite  bi-axial  strain  states;  lattice  mismatch  is  -1.1%  on  LAO  and  +1.8%  on  STO.  To  probe 
the  response  of  the  Ni  3 d  orbital  occupations,  we  employed  x-ray  absorption  spectroscopy  (XAS) 
with  linearly  polarized  lights,  i.e.  x-ray  linear  dichroism  (XLD),  at  the  Ni  L-edge.  Figure  8  shows 
the  normalized  absorption  spectra  taken  by  the  fluorescence  yield  mode  at  room  temperature  for 
in-plane  (H=L  _L  c)  and  out-of-plane  (V=E  II  c)  polarizations  under  the  two  opposite  strain  states. 
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Figure  8:  Normalized  absorption  spectra  at  the  NiL-edge  taken  in  the  bulk-sensitive  fluorescence  mode  with 
varying  photon  polarization  (see  legend)  for  films  under  compressive  (red)  (a)  and  tensile  (blue)  (B)  strain,  on  LAO 
and  STO  substrates,  respectively.  The  black  curves  in  (c)  represent  the  bulk  (ceramic)  LNO  data.  (RIGHT) 
Schematic  orbital  level  diagram  for  bulk  LNO  and  the  anticipated  strain-induced  orbital  polarization  effect  on  the  eg 
doublet. 

As  seen,  each  polarization  averaged  spectra  is  representative  of  octahedrally  coordinated  Np+ 
and  in  good  agreement  with  the  spectrum  taken  on  bulk  LNO  (see  Fig.8(c)).  There  are  distinct 
differences,  however,  in  the  polarization  dependence  between  films  under  compressive  and 
tensile  strain.  Specifically,  as  showed  in  Fig.  8(a),  the  absorption  for  the  out-of-plane  polarization 
is  shifted  by  -100  meV  lower  in  energy  than  the  in-plane  polarization  for  LNO  under 
compressive  strain  on  LAO.  This  dichroism  indicates  a  small  eg  band  splitting  between  the  Ni- 
2  2  2  2 

derived  x  -  y  and  3 z  -  r  orbitals  due  to  the  distortion  of  the  spherical  charge  distribution 
around  the  Ni  sites  by  non-uniform  Ni-O  bond  lengths.  The  sign  of  the  XLD  implies  that  the  in¬ 
plane  Ni-O  bonds  are  moderately  compressed  relative  to  the  out-of-plane  apical  Ni-O  bonds 
along  the  c-axis  in  agreement  with  the  model  of  SIOP  concept. 
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Figure  9:  Different  types  of  distortions  for  accommodating  lattice  mismatch.  Dashed  diamonds  stand  for  undistorted 
octahedral  network,  while  the  solid  ones  represent  the  distorted  network.  Arrows  denote  the  distortion  motions,  a) 
compressed  octahedron  in  accordance  with  SIOP.  b)  cooperative  octahedral  rotation,  c)  breathing  mode. 


For  tensile  strain  with  the  same  substrate-film  orientation,  the  SIOP  model  predicts  an  inversion 
of  the  orbital  polarization  and  the  population  of  the  dx2_y2  orbital;  however,  the  polarization 

dependence  of  our  films  on  STO  is  remarkably  different.  As  see,  the  orbital  dichroism  is 
practically  absent.  This  result  implies  that  the  uniform  Ni-O  bond  lengths  and  the  highly 
symmetric  crystal  field  remain  seemingly  intact  despite  the  large  (~  +2%)  tensile  strain. 

These  findings  are  surprising  since,  according  to  the  SIOP  framework,  the  larger  magnitude  of 
the  strain  from  STO  than  LAO  should  also  produce  a  large  orbital  polarization  albeit  with  an 
opposite  sign.  The  violation  of  the  SIOP  model  under  the  tensile  strain  demonstrates  an 
asymmetrical  orbital  response  to  the  bi-axial  epitaxial  strain,  and  is  indicative  of  other  efficient 
competing  mechanisms  in  accommodating  the  lattice  mismatch.  One  possibility  is  the  octahedral 
rotation  which  can  modulate  the  pseudocubic  lattice  parameter  by  changing  the  Ni-O-Ni  bond 
angle  instead  of  the  Ni-0  bond  length.  Density  functional  calculations  within  the  local-spin 
density  approximation  (LSDA)  plus  Hubbard  U  method  were  also  computed  based  on  our 
experimentally  determined  lattice  parameters.  The  resulting  structural  optimizations  are  indeed 
consistent  with  our  findings  by  XLD.  While  the  structure  on  LAO  consists  of  a  single  Ni  site 
with  different  in-plane  and  out-of-plane  Ni-O  bond  lengths  that  lift  the  formerly  bulk-like  cubic 
Ni-site  symmetry  and  the  orbital  degeneracy,  it  show  a  pronounced  pattern  of  octahedral 
rotations  on  STO.  Moreover,  the  tensile  strain  also  stabilize  an  additional  intriguing  breathing 
mode  imposed  on  the  Ni06  octahedral  network  to  accommodate  the  lattice  mismatch.  This 
‘breathing’  octahedral  distortion  induces  a  Ni-O  bond  length  disproportionation  by  splitting  the 
equivalent  Ni  sites  into  Ni(l)  and  Ni(2)  that  are  located  at  the  center  of  nearly  ideal  NiOg 

octahedra  with  different  averageNi-O  bonds  lengths  of  2. 00 A  andl  .93 A  respectively. 
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